Spectroscopy of atomic rubidium at 500 bar buffer gas pressure: approaching the 
thermal equilibrium of dressed atom-light states 
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We have recorded fluorescence spectra of the atomic rubidium D-lines in the presence of several 
hundreds of bars buffer gas pressure. With additional saturation broadening a spectral linewidth 
comparable to the thermal energy of the atoms in the heated gas cell is achieved. An intensity- 
dependent blue asymmetry of the spectra is observed, which becomes increasingly pronounced when 
extrapolating to infinitely high light intensity. We interpret our results as evidence for the dressed 
(coupled atom-light) states to approach thermal equilibrium. 

PACS numbers: 32.80.-t, 05.30.Jp, 32.70.Jz, 42.50.Fx 



In usual optical spectroscopy experiments, light 
sources along with spontaneous decay drive the inves- 
tigated sample far from thermal equilibrium This is 
most evident in room temperature atomic physics exper- 
iments, where the observed spectral linewidths in energy 
units are many orders of magnitude below the thermal 
energy k B T (~ 1/40 eV for T = 300 K). Thus, apart 
from a comparatively small and symmetric broadening 
arising from the known Doppler effect, the excitation pro- 
file is largely independent of the statistical distribution 
function. In cold atom experiments, spectral linewidth 
and thermal energy can be comparable. However, the 
lack of a sufficiently fast thermalization process has so 
far here prevented thermal equilibrium of coupled atom- 
light states to be a useful concept. 

We here report on experiments, in which up to 500 bar 
of argon and helium respectively buffer gas pressure in- 
duce a linewidth of a few nanometers for the rubidium 
D-lines. At high optical power of the exciting continuous- 
wave laser source, the fluorescence spectra are broadened 
by additional power broadening to values exceeding the 
thermal energy ksT in the heated gas cell. In this regime, 
we observe a strong blue asymmetry of the spectra. The 
spectral asymmetry increases further when extrapolating 
our data towards infinitely high excitation intensity. We 
interpret our results as evidence for the coupled atom- 
light states (" dressed states" ) to approach thermal equi- 
librium, with the thermalization being due to frequent 
rubidium- buffer gas collisions. Thermal equilibrium is 
e.g. a prerequisite for possible BEC-like phase transitions 
of coupled atom-light quasiparticles (polaritons) Q. 

Our experiment benefits of collisional aided excitation, 
which allows to saturate the atomic transition at large 
detuning with a continuous-wave laser source. Pressure 
broadening of atomic spectral lines is a long investigated 
topic [1, U H, @, 0> H, M, H3] • For large laser detunings 
(S> 1/tcoU, where t co u denotes the collisonal duration), 
the impact limit ceases to be fulfilled, and atomic lines of- 
ten become asymmetric. A theoretical description of the 
absorption probability requires knowledge of the molec- 



ular potential curves of the collisional partners llj, and 



can be expressed e.g. in terms of modulated dipole ap- 
proaches [12| • The effect of collisional redistribution has 
shown the influence of state changing collisions In 
an interesting experiment, deviations from the Einstein 
coefficients for absorption and stimulated emission have 
been observed for far off-resonant excitation in a colli- 
sionally broadened system 14[ . Dressed state approaches 
have proven to be an elegant way to treat collisional pro- 
cesses in the presence of laser radiation [l5| . 

For a theoretical description of the transition to ther- 
mal equilibrium of dressed states, we give a simple model 
based on thermodynamic arguments. Consider a two- 
level system with ground state \g) and excited state |e) 
coupled to a laser field of frequency u. When tuned into 
resonance, the laser field connects the states \g,n + 1) 
and |e, n) respectively, where the first quantum number 
denotes the internal atomic state and the second one the 
photon number. To simplify the analysis, let us restrict 
ourselves to the case of large detunings 8 = uj - ui atom , 
with \b\ ^> ilRabi, so that in the absence of collisions both 
the state mixing and the AC Stark shift can be neglected. 
The energy splitting between dressed states is thus sim- 
ply S, as shown in Fig. la. On the other hand, when a 
buffer gas atom approaches the energy levels are shifted, 
and efficient excitation and a transfer between dressed 
state energy levels is provided when the laser field tran- 
siently becomes resonant during the collision, yielding an 
example of collisional aided excitation. Let c 9 _> e (<J) and 
c e ^ g (5) denote the rate constants for collisional transfer 
between the dressed state levels, which for \5\ Sljiabi 
can be shown to be proportional to the laser intensity 
[jjj . For a finite spontaneous decay rate T, also a cou- 
pling between states \e, n) and \g, n) of different dressed 
state manifolds occurs. 

We are interested in the properties of such a two level 
system subject to laser radiation and coupled to a ther- 
mal bath. If we assume our system to be optically thin, 
we can trace over the photonic quantum numbers, and 
obtain the following Boltzmann-like rate equations 
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FIG. 1: (a) Scheme of levels and rates for a two-level dressed 
state system approaching thermal equilibrium in the far de- 
tuned limit, (b) Upper state population as a function of laser 
frequency, where a Lorentzian lineshape was taken for the 
detuning dependent excitation rate c(5): For c(fcsT) <C T, 
i.e. small linewidth, which yields a lineshape resembling that 
of c(S) (dotted line), c(k B T) ~ T (dashed line) and c(5) > T, 
i.e. thermal equilibrium of dressed states (solid line). 

f e = Cg^ e (S)fg - C e ^g(6)f e - Tf e (la) 

fg = C e ^g(S)f e - Cg^ e (S)fg + Tf e , (lb) 

where f e and f g are the corresponding populations of 
states |e) and \g). To derive a relation between the 
rates c g —, e and c e _» g , let us consider the equilibrium so- 
lution (/ e = f g = 0) in the absence of spontaneous decay 
(i.e. r — > 0). In this limit, we expect that usual Boltz- 
mann statistics can be applied to this two-level system, 
which yields f e = e & ' ksT / (l + e 6 / kBT ) = l/(l+e- 5 / kBT ), 
which resembles the Fermi-Dirac distribution function 
due to only two available states. From Eqs. 1 we obtain 
f e = 1 — f g = 1/(1 + c e —>g/cg^> e ) so that the following 
relation of the rate constants must apply: c e —, g /c g —> e — 
e -5/k B T ^ ag a i g0 s ] 10wn m pjj] Note that the sign of the 

laser detuning in this formula depends on the definition 
of 6. Physically, the increased rate towards a popula- 
tion of the energetically lower dressed state is due to the 
larger phase space that becomes available to the reservoir 
with the associated energy gain of \8\. The result is also 
obtained when considering the variation of the available 
density of states in Fermi's golden rule. 

Spontaneous decay drives the dressed state system out 
of thermal equilibrium. The full Eq. 1 with f e = f g = 
yields the following stationary solution in our model: 
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The condition for this formula to reduce to the usual 
Fermi-Dirac distribution of a two-level system at large 
transition linewidth and laser intensity is that c(5) 3> T 
is fulfilled at all relevant detunings 8. On the other hand, 
for k B T > \5\ we obtain f e ~ (c(S)/T)/(l + 2c(S)/T), 
which corresponds to the usual saturation dependence of 
a two-level system 
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Fig. lb shows a comparison of 
line profiles in the different regimes, where a Lorentzian 
excitation profile was assumed for the detuning depen- 



dent excitation rate c(5) . The dotted line gives the spec- 
trum in the limit of a linewidth of c(5) considerably be- 
low k B T, resulting in a lineshape resembling that of c(<5), 
while the solid line gives the case of a linewidth con- 
siderably above the thermal energy (solid line). Also, 
an example for the intermediate regime of comparable 
linewidth and thermal energy is given (dashed line). A 
theoretical prediction of the detuning dependent excita- 
tion rate c(8) in the here investigated dense buffer gas 
regime is a highly nontrivial task, especially since we 
expect multiple particle collisions to play an important 
role. On the other hand, if full thermal equilibium can 
be achieved for all used laser detunings, the predicted 
lineshape reduces to a Fermi-Dirac profile independently 
of the form of c(5). 
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FIG. 2: Scheme of the experimental setup. 

A scheme of our experimental setup used to investi- 
gate rubidium atoms at high buffer gas pressures and 
intense laser radiation is shown in Fig. 2. A rubidium 
filled stainless steel oven with volume of a few cm 3 and 
wall thickness of roughly a cm is used as a pressure cell. 
Optical radiation is coupled into the chamber through 
sapphire windows. Buffer gas is filled into the chamber 
through an attached valve with pressures up to 230 bar 
(at room temperature). Due to technical reasons, the 
pressure is measured outside the chamber. During op- 
eration, the cell is heated to 260 °C, which results in 
a vapour pressure limited rubidium number density of 
~ 1.0 • 10 16 /cm 3 [ill]. In the sealed chamber, the buffer 
gas pressure then increases up to approximately 400 bar 
for helium and 500 bar for argon. 

Tunable laser radiation for excitation of the atoms is 
provided by a Ti:sapphire ring laser. Laser frequency 
scanning over the broad atomic absorption profile is done 
manually within the operation range of the used laser 
mirror set of 350-395 THz. The laser frequency is mea- 
sured with a wavemeter. With a Fraunhofer achromat, 
the laser radiation is focused into the rubidium chamber 
to a beam radius of ~ 3 /im, which at a maximum op- 
tical power of 300 mW results in an intensity of 1.1T0 8 
mW/cm 2 . The beam focal plane within the cell is placed 
near the sapphire window. To selectively record the 
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atomic fluorescence only from regions of high laser in- 
tensity, we furthermore spatially filter both the incident 
optical beam and the outgoing fluorescence with pinholes 
in a confocal geometry. The recorded fluorescence power 
values exhibited short-term variations that became in- 
creasingly severe for a focal spot very close to the cell 
window, which attribute to thermal fluctuations. For our 
experimental spectra, the mean fluorescence power value 
averaged over 30 seconds is taken for every data point. 

In initial experiments, we have measured the lifetime 
of the rubidium 5P state at high argon buffer gas pres- 
sures. For this measurement, an acoustooptic modulator 
(not shown in Fig. 2) was included in the beam path, 
which allowed for a rapid shutting of the exciting optical 
radiation for up to ~ 30 mW beam power. The sub- 
sequent decay of the atomic fluorescence was monitored 
with a photomultiplier tube. At low buffer gas pressures 
(< 10 bar) we observed fluorescence decay times up to 1 
/is, which we attribute to repopulation from higher states 
due to energy pooling [17| . More direct evidence for en- 
ergy pooling is obtained from the clearly visible blue flu- 
orescence at those pressure values. At higher buffer gas 
pressures the fluorescence signal decayed with a time con- 
stant of 40 ns. If we account for the time constants of 
our acoustooptic modulator and the photomultiplier de- 
tection, this value is in good agreement with the lifetime 
value of about 27 ns obtained in earlier measurements 
[l8j |. We conclude that energy pooling to higher state 
excitation is suppressed at sufficient buffer gas pressure. 
Further, the lifetime of the 5Pi/2 and 5P3/2 states are, 
within our measurement accuracy of 5 ns, unchanged 
to argon pressures of 200 bar, which is the maximum 
achieved pressure value in these early data sets. At a col- 
lisional rate of 10 11 /s, the atoms experience more than 
10 3 collisions within a natural lifetime. This confirms the 
remarkable elasticity of excited atoms with rare gas colli- 
sions [l9| . and is a prerequisite for the thermalization of 
dressed states in our experimental scheme. Note that the 
excited state lifetime here is some four orders of magni- 
tude above the ~ps lifetimes observed in semiconductor 
exciton systems |2j. 

Typical fluorescence spectra recorded under high pres- 
sure buffer gas conditions for variable light intensities are 
shown in Fig. 3a (argon buffer gas) and Fig. 3b (helium 
buffer gas). While the observed lineshapes somewhat 
differ for both used buffer gases, indicating the 'tran- 
siently chemical' nature of optical collisions in a regime 
beyond the impact limit, let us draw our attention to 
some general features for both spectra. At moderate 
optical power (P ~ 25 mW) the linewidths are clearly 
below the thermal energy k B T (~ 11 THz at T = 530 
K). At higher power levels the linewidths increase to val- 
ues around fcgT. Noticeably, the red side of the Dl- 
line tends to saturate to a comparatively low fluorescence 
level, while at the blue side of the D2-linc much higher 
fluorescence is observed, though in the far wings even at 
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FIG. 3: Fluorescence spectrum of the rubidium D-lines at 500 
bar argon (a) and 400 bar helium (b) buffer gas pressure for 
different optical beam powers. The small drawings on the top 
indicate the population of the dressed states on the red and 
blue side of the electronic transition respectively. 



maximum power levels (300 mW) the transition is less 
saturated. We attribute this to the dressed state system 
approaching thermal equilibrium for high power levels. 
The resonant Rabi coupling at full optical power is of 
order 100 GHz, so that for detunings of order fc^T the 
mixing of ground and excited states is still small. For 
the sake of simplicity, let us restrict the discussion to 
the red side of the Dl-line and the blue side of the D2- 
line, as in these limits the influence of the upper state 
fine structure splitting is smallest. While on the red side 
of a two-level system the state \g, n + 1) is energetically 
below |e,n), leading to an enhanced ground state pop- 
ulation in thermal equilibrium, on the blue side of the 
transition the dressed excited state component is lowest 
and expected to be favored in equilibrium, as indicated in 
the boxes in the top of Fig. 3. At low optical power lev- 
els, the upper state spontaneous decay drives the coupled 
atom light system away from equilibrium. For a more de- 
tailed analysis, Fig. 4a shows the observed fluorescence 
as a function of the optical power for three different laser 
frequencies for the data with helium buffer gas. These 
saturation curves have been fitted with the theoretical 
prediction of Eq. ([2]) , where a constant prefactor and the 
ratio of c(S) and the laser intensity were left as free fit 
parameters. This fitting procedure assumes that the de- 
tuning dependent excitation rate c(S) is linearly depen- 
dent on the laser intensity, an assumption in agreement 
with present theories of optical collisions, whose valid- 
ity in the yet largely unexplored high pressure buffer gas 
regime with large multiparticle collisional rate however 



remains to be tested [21(. Near a line center (crosses), 
the curve saturates already at relatively low power level 
to an intermediate fluorescence level. For significant de- 
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FIG. 4: (a) fluorescence of rubidium atoms at 400 bar helium 
buffer e SlS £LS cL function of optical power for different laser fre- 
quencies. The experimental data has been fitted with the ex- 
pected saturation dependence of Eq. @. (b) Expected atomic 
fluorescence at infinite laser power, obtained from extrapola- 
tion curves as in (a), versus laser frequency. The shown error 
bars were obtained from the extrapolation fits. The data has 
been fitted with a Fermi-Dirac distribution, assuming thermal 
equilibrium of the two-level dressed state system. 



tuning the optical power at which saturation is achieved 
is higher, as visible both for blue (dots) and red (circles) 
detunings. Interestingly, the saturation level of fluores- 
cence is clearly different in both cases, with the blue (red) 
detuning leading to largest (smallest) values. 

Fig. 4b shows the value of fluorescence extrapolated 
to infinite beam power, as taken from such curves, as 
a function of laser frequency. This is to provide an es- 
timate for the upper state population at thermal equi- 
librium. The data has been fitted with a Fermi-Dirac 
distribution f e (5) = 1/(1 + e ~ <5 / fcsT ) expected for the 
upper state population of the two level dressed state sys- 
tem, where the detuning 5 was measured relatively to 
the center of the D-lines. The agreement between the- 
ory and experiment is quite reasonable, which supports 
the assumption of the dressed states to approach ther- 
mal equilibrium at high laser power. An issue remaining 
to be clarified is that a best fit is obtained with a tem- 
perature T = 162(58)K, which is significantly below the 
cell temperature. For a comparable fit with the argon 
buffer gas spectral data, a temperature of 270(91)K was 
obtained. At present, the origin of this disagreement of 
temperatures is not resolved. A more detailed rate equa- 
tion model of partial thermal equilibrium should include 
the rubidium fine structure and spatially inhomogeneous 
saturation. On the other hand, an interesting question is 
whether laser cooling (heating) processes of the thermal 
gas here occur on the red (blue) side of the transition due 
to the energy loss (gain) during collisional aided excita- 
tion, as was proposed in [20| . If the spontaneously emit- 
ted photons are energy shifted by an amount comparable 
to fcflT respectively to the incident photon wavelength, 
cooling may be possible in such pressure broadened sys- 
tems, comparable to results obtained on the optical cool- 



ing of solids [22J . 

To conclude, we have recorded spectra of rubidium 
atoms in a collisionally broadened regime interpolating 
between usual atomic physics gas phase and solid/liquid 
phase conditions. With additional saturation broaden- 
ing, the spectral linewidth approaches the thermal en- 
ergy. The saturation dependence of the line profile is 
interpreted as evidence for the onset of thermal equilib- 
rium of dressed atom-light states. 

For the future, it would be interesting to study pos- 
sible novel laser cooling mechanisms of high pressure 
atomic gases based on collisionally aided excitation. A 
quite different perspective could include a possible BEC- 
like phase transition to a condensed atom-light polariton 
phase. For our present apparatus, this would require an 
increase of the optical intensity and an optically thick 
atomic ensemble, as could be realized with a cavity- or 
optical waveguide based system. 

We acknowledge helpful discussions with N. Schopohl 
and R. Walser. This work has been funded within a col- 
laborative transregional research center (TR 21) of the 
DFG. 
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